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O |. Abstract 

We consider a unified model, called ancestor model, associated with twisted trigonometric R matrix which model 
leads to several descendant integrable lattice models related to the \J q i/2(sl 2 ) symmetry. Boundary operators 
compatible with integrability are introduced to this model. Reflection and dual reflection equations to ensure 
integrability of the system are shown to be same as the untwisted case. It follows that underlying symmetry of the 
ancestor model with integrable boundaries is identified with the q-deformed analogue of Onsager's symmetry. The 
transfer matrix and its related mutually commuting quantities are expressed in terms of an abelian subalgebra in 
the g-Onsager algebra. It is illustrated that the generalized McCoy- Wu model with general open boundaries enjoys 
t-H , this symmetry. 
OV 
O 
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1 Introduction 



Unifying scheme to rich varieties of models has played important roles to develope our concepts and perspectives 
in all scientific fields. Especially in theoretical physics and mathematics, most of successes of such attempts is 
deeply connected with findings of underlying symmetries, which brings beauty and simplicity into relief. Although 
such symmetries may be hidden, in many cases they are codified by groups which are generated by Lie algebra and 
its generalizations including the Kac-Moody algebra and quantum algebra. One of the most direct application in 
this subject is in the theory of nontrivial completely integrable models (continuum or lattice) in 1+1 space time 
dimensions. Their models have an (in)finite set of sufficient mutually commuting conserved quantities that all 
physically relevant information can be extracted without any approximation in principle. Therefore the interesting 
features of integrabilities have made profound influence between many branches of physics and mathematics, for 
example (super)string theory and condensed matter systems in physics and quantum Lie algebra, g-orthogonal 
polynomials, g-difference equations and knot theory in mathematics [1, 2, 3, 4, 5]. 

Since the seminal works by Cherednik [6] and Sklyanin [7] , quantum boundary integrable models have received a 
lot of attention in the context of systematic studies on the XXZ anisotropic open spin chain models [8, 9, 10, 11, 12], 
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two dimensional boundary integrable quantum field theory [13, 14, 15, 16, 17, 18, 19, 20], the open string/spin chain 
sector of AdS/CFT correspondence in super string theory [21, 22, 23, 24] and so on. In recent years, based on the 
Sklyanin's dressing method, a breakthrough emerged [25, 26] from the study of the algebraic structure encoded in the 
reflection equation for certain boundary conditions and U^i^s^) trigonometric i?-matrix. Furthermore successive 
studies in this line revealed that the transfer matrix of the model (including the XXZ open spin chain model with 
general boundary conditions) is generated by an abelian subalgebra of the q-Onsager algebra [27, 28]. The g-Onsager 
algebra is also referred as the tridiagonal algebra, which is first introduced by Terwilligcr in the context of the P and 
Q- association schemes related to the Askey scheme of orthogonal polynomials in the mathematical side [29]. The q- 
Onsager algebra can be viewed as an extension of the q-Serre relations and a g-anlogue of the Dolan-Grady relations 
[30] that define the Onsager algebra [31]. It is worth emphasizing that the Onsager algebra not only have played 
crucial roles in planar Ising model, XY model and superintegrable chiral Potts model [31, 32, 33, 34, 35, 36, 37], but 
so much of what we understand in generality was contained in the Onsager's solution that led us to the discovery 
of the Yang-Baxter equation (the star triangular relation) [38]. It seems that the emergenece of the g-Onsager 
algebra opens the new possibility of studying massive continum or lattice quantum integrable models. 

New possibilities to deform quantum integrable (spin chain) models using the Drinfeld twist [39] have already 
been considered [44, 45]. We employ their ideas to construct much wider class of quantum integrable lattice systems 
and introduce integrable boundaries for the ancestor model [44] related to the V q (sl2) case. Then matrix elements 
of the associated quantum Lax operator is made of t (or #)-deformed extension of the trigonometric Sklyanin 
algebra [44], which describe several quantum integrable models related to the twisted U g (s?2) in a unified way. 
The main purpurse in the article is to identify underlying symmetry in the ancestor model with open integrable 
boundary conditions, by applying the Sklyanin's dressing method with the Lax operator. 

This article is organized in the following way: in the section 2, a brief formulation of quantum inverse scattering 
method (QISM) for quantum integrable systems and its integrable boundary extension is summarized. The Lax 
operator for the twisted V q (sl2) case is uniquely determined except for the gauge degreee of freedom up to the 
Laurant polynomials of degree 1 in the section 3. To utilize the Sklyanin's dressing method, we need to find 
the inverse representation of the Lax operator. To find the inverse operator can be performed by introducing 
a new operator. It is shown that a new operator unentangles nontriviality of the twisted transformation to the 
untwisted algebra. In the section 4, it is shown that reflection algebra resulting from the reflection equation for the 
twisted U g (s?2) case is the same as the untwisted one (the extended trigonometric Sklyanin algebra). As shown in 
Baseilhac and Shigechi [40], it is clear that the g-Onsager symmetry appears in the reflection algebra. We explicitly 
construct the Sklyanin's dressed solution with the c-number solution by Goshal and Zamoldochikov. As expected, 
the foundamcntal operators for N = 1 dressed solution are satisfied with the Askey- Wsilson relations. Through 
the same procedure in the articles [27, 28], general N dressed solution is explicitly constructed. The fundamental 
operators are satisfied with the tridiagonal relation, so that the transfer matrix enjoys the g-Onsager symmetry. As 
described by linear combinations of generators of the q-Onsager symmetry, the transfer matrix is made from the 
abelian subalgebra of the q-Onsager algebra. As the application of this construction, we derived the Hamiltonian 
of a generalized McCoy- Wu model [41, 44] with general boundary conditions in the section 5. The model may be 
considered as a generalization of XXZ spin chain with Dzyyaloshinsky-Moriya interactions [42] . All of the integral 
motions arc identified as an abelian subalgebra of the q-Onsager algebra. The discussion is given in the section 6. 

2 Formulation on quantum integrable lattice systems 

This section serves as a brief introduction by quantum inverse scattering method [43] for quantum lattice systems 
with periodic and open boundaries compatible with integrability. 

Complete integrablity of quantum systems takes its roots on the existence of (in)finite set of mutually con- 
served quantities to be involutive. Within quantum inverse scattering method, this is encoded in the existence of 
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monodromy matrix T(it) satisfying the quantum Yang-Baxter equation(QYBE) 

R\2{u,v) T (u) T (v) =T (v) T (u)R 12 (u,v), (1) 

1 2 

where T= T ® J, T= I <S> T, which equation is in auxiliary spaces V\ ® V 2 and a full Hilbert space % acted by the 
monodromy matrix T. Taking the trace of the above equation, we obtain 

Tr 12 (T («) T («)) = Tr 12 (R^(u,v) T («) T (u)i? 12 (u, «)) = Tr 12 (T («) T (u)), (2) 

12 12 

where Tri 2 (T (u) T (v)) defines the trace over V\ ® V 2 . By using the definitive relation Tri 2 (T' (11) T (^)) — 
Tr(T(tt))Tr(T(u)), it is shown that the transfer matrix t(u) = Tr(T(u)) is satisfied with the trivial commutation 
relation 

[t(u),t(v)] = 0, (for any variables u and v). (3) 

The above equation guarantees that the transfer matrix posecscs mutually commuting quantities /2/c+i as the from 
t( u ) = 5Z-? 7 2fe+i('u)/2fe+i- As the result, the integrability of the quantum systems is ensured in the quantum version 
of the Liouville sense defined in the context of the classical case. 

For TV-site periodic lattice models, the monodromy matrix T(u) is discribed by 

T(u) = L N {u)L N ^{u) ■ ■ ■ L 2 (u)Li(u) (4) 

in terms of the Lax operator Li (u) at the i-th lattice site acting on a local Hilbert space Hi . Requiring ultralocality 

1 2 

for the Lax operator, [Li, Lj] — (i ^ j), the QYBE leads to 

Ri2(u,v) Li (u) Li (v) =L (v) L (u)R 12 (u,v), (i = l,---,N) (5) 

12 3 

as a consequence of its local version. Supposed that the associativity for the triple product Li (u) Li (v) Li (w) 
holds, the equation in effect turns into an independent equation for the matrix R(u) known as the Yang-Baxter 
equation (YBE) 

R 12 (u,v)R 13 (u,w)R 23 (v,w) = R 23 (v,w)R 13 (u, w)Ri 2 (u,v). (6) 

For ultralocal lattice models that we are intrested in, the matrix R\ 2 (u, v) takes the form of R^uv^ 1 ), i.e. the 
matrix elements of R are given through functions of the spectral parameters uv^ 1 . 

For quantum integrable systems with open boundaries, each set of boundary conditions is associated with a 
choice of boundary operators K± (u) . Without spoiling the nice algebraic structure and analytic properties of bulk 
integrability, the adequete conditions for the boundary operators K± (u) compatiable with integrability, respectively, 
are obtained from the following equations known as the reflection and dual reflection equations [6, 7, 46]: 

^(uw -1 )^-^)®!)^* 3 ^)^® - {l®K-{v))R 12 {uv){K-{u)®l)R t J 2 t2 {uv- 1 ), (7) 

R l2 {u- 1 v){K i ^{u) ® \){M~ l ® 1)-R* 2 * 2 (q~ 1 u~ 1 v~ 1 )(M ® 1)(1 ® K^(v)) 

= {\®K t2 {v)){M <g> l)i? 12 (q- 1 u- 1 w- 1 )(M- 1 ® l)(K*2{u) ® l)^* 2 ^"^), (8) 

where the symbol U denotes transposition in the i-th auxiliary space. Furthermore the matrix M is determined by 
the following relation [47] 

{{{^(w)}" 1 }* 2 }" 1 = ® M )Ri2{qu){l ® M)- 1 , M* = M, (9) 

where ({u)l[ = i?i 2 (u)i? 2 i( u_1 ) with the unit matrix I. Once finding out a solution of the reflection equation (7) 
for K_(u), one can verify that K + (u) matrix defined by 

K + (u) = Kliq- 1 / 2 ^ 1 )]^ (10) 
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satisfies the dual reflection equation (8). The transfer matrix for the TV-site lattice model with the integrable open 
boundary conditions is defined by 

t(u) = Tr(K + (u)K^_ N \u))), (11) 

where 

K {N) (u) = ((L N (uv N ) ■ ■ ■ Li(««i)) K_{u) (L^ 1 (u- 1 v 1 ) ■ ■ -L^u" V))) • (12) 

Interestingly the dressed boundary operator (u) is also a solution of the reflection equation. This permits us 
to make a sequence of the operator- valued solutions from one solution of the reflection equation. This construction 
is known as the Sklyanin's dressing method. Thanks to these algebraic relations, the transfer matrix commutes for 
any spectral parameters u and v. 

[t(u),t(v)]=0, (13) 

which is enough to ensure the integrabiity of the system even with open boundaries. 

To construct much wider class of integrable models starting from a solution of YBE for the i?-matrix, let's look 
into a transformation called twist introduced in the article [39]. Provided that the twisting operator T satisfies the 
conditions 

i?12(u)^ 7 12^ : 23 = ^23^1zR\2(,u) , -^l 2 -^13 -^2 3 = FiaFlzFll, Tij—T^ , (14) 

we obtain the twisted i?-matrix 

Ri2(v) = T^R 12 {u)F^ (15) 

from the original one R\ 2 (u). Corresponding to the transformation, twisted Lax operator and twisted boundary 
operators are obtained. 



3 Lax operator related to twisted trigonometric i?-matrix 



The Yang-Baxter equation (6) restricts the solution of the i?-matrix to integrable models. The most simplest 
solution of YBE acting on the two dimensional auxiliary spaces V\ <8> V2 is the trigonometric solution related to 
\] q i/2{sl 2 ) symmetry: 

( a(u) \ 
b(u) c 
c b(u) 



R 12 (u) 



a(u) 



V 



(16) 



/ 



where a(u) = q 1 / 2 u — q 1 I 2 ^ 



bin) 



c = q 1 ! 2 — q 1 / 2 . Starting from the trigonometric solution 



(16) of YBE, one can produce the twisted trigonometric i?-matrix in terms of a suitable representation T12 
e i6(rT3/2®z-z®a 3 /2) ag twisted operator, where Z is the central charge. As a result, we obtain 



R 12 {u) = 



( a ( u ) 



V 



tb(u) c 
c t- x b(u) 



a(u) 



(17) 



where t = 



n -2i0Z 



The twisting procedure enables the Lax operator to take the different twist parameter t n at 



each lattice site. Then the related twisted Lax operator L n {u) at the n-th lattice site takes 2x2 matrix form in 
the auxialiary space as follows: 



L n (u) = 



UT-i 



721 



Tl2 



(18) 
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These operators Tf 1 and (i,j — 1,2) satisfy the i-deformation of the extended trigonometric Sklyanin algebra 
[44]: 

t±t 12 = i-V^Wf, r^ris = t- 1 « =Fl/2 ri2T 2 ± , 

T 1 ± T 2 1 = t„g Tl/2 T 2 lT 1 ± , T^T12 = t n q ±1/2 T 2 lT^ 7 

t»T 2 lTl2 - t~ 1 T 12 T 2 l = c(t^T2 - T^T^). (19) 

The coproduct structure for this algebra is found by the elegant formulation by Faddeev-Reshetikhin-Takhtajan 
[48]. Symmetric and nonsymmetric realizations of the algebra can generate several descendant lattice models 
without limiting procedures, for examples the generalized McCoy- Wu model and ^-deformation of such models 
as the Liouville lattice model, q-oscillator model related to the sl q (2) and sZ q (l, 1) [49], the derivative non-linear 
Schrodinger equation, lattice version of the massive Thirring model, lattice sine-Gordon and their hybrid models 
[50] and so on. Although one could explicitly construct the transfer matrix for these integrable models with periodic 
boundary conditions, the derivation is out of our intrests in the present aricle. 

Our main aim in this article is to introduce intgrable open boundaries for these models and construct the 
transfer matrix and mutually commuting quantities in terms of generators in the q-Onsager symmetry. For our 
purpurse, it is useful to factorize the above opearators (19) in terms of new operators 

T t = TtTg,Ti 2 = t~ 1/2 f 12 T g , T21 = # 2 f 2 lT fl , (20) 

with the following algebraic relations 

Tgf 12 =t- 1 f 12 Tg, Tgf 2 l=t n f 2 lTg, [<T g , T f } = . (21) 

Plugging the factorization of the operators into Eq.(19), they reduce to the algebraic relations on the non-twisted 
extended trigonometric Sklyanin algebra: 

nfi2 = <7 ±1/2 fi 2 f 1 ± , f 2 fi 2 = q Tl/2 f 12 f 2 , 

T\f 2 \ = <7 Tl/2 f 2 if 1 ± , f 2 f 2 i = (7 ±1/2 f 2 if 2 ± , 

[751, fi 2 ] = c(f+f 2 - - ff r+), (22) 

It turns out that the introduction of these new operators unpicks the entanglement of twisted quadratic algebra. 
As usual, it is shown that there exist five Casimir operators w±, wo x , wo 2 , w: 

f±f± = w ±, frff = w 0t , (i = 1, 2), 
712721 - 9 1/2 fff 2 + - q- 1/2 fff 2 = f 2 ifi2 - q~ 1/2 f^f^ - q 1/2 fff 2 = w, (23) 

with the relation w-w + = w 0l wo 2 - 

The Lax operator (18) is re-expressed in terms of these algebraic elements: 

T I \ { U7f +M _1 f 1 + t^ 1/2 fl 2 \ , . 

Ln{u) = 1/2 . T g . (24) 

Y trl T21 UT 2 + U 1 T 2 r J 

As these consequences, we obtain the operator L n (u) 
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which is proportional to the inverse of the Lax operator. In fact, it is easy to check 

L n (u)L n (u) = p(u)I, p(u) = w — (g _1 / 2 w_« 2 + q^-^w+u^ 2 ) 
with the algebraic relations (22) and the Casimir operators (23). 



(26) 



4 Reflection equation and dressing K matrix related to twisted R- 
matrix 

In the previous section we derived the fundamental parts to construct dressed if-matrix, although solutions for 
the reflection equation and the dual reflection equation are still left. We explicitly construct thier solutions in this 
section. 

Before finding solutions of the reflection equation and dual reflection equation for boundary operators K±(u), 
it is worthwhile writing out all elements of the reflection equation. Supposed that 

/ A(u) B{u) \ 
A " W - ^ C(u) D(u) ) ' 

the reflection equation for K_ reduces to the sixteen algebraic equations: 

(1) a_c(BC' - B'C) + a_a + [A,A'} =0, 

(2) a_c(CB' - C'B) + a_a+[L>, D'} = 0, 

(3) b-b+[A, D'} + £ 2 [D, D'} + ca+{CB' - C'B) = 0, 

(4) b-b + [D,A l ] + ~c 2 [A,A']+£a + {BC l - B'C) = 0, 

(5) cb+(DA' - D'A) + b^c{AA' - DD') + b-a+[B, C] = 0, 

(6) cb+(AD' - A'D) + b^c{DD' - AA') + fe_a+[C, B'] = 0, 

(7) b-b+AC + c?DC + ca+CA' - a^a+C'A - a-cD'C = 

(8) b-b+DB' + c 2 AB' + m+BD' - ci-a+B'D - a-cA' 'B = 0, 

(9) b-b+B'A + c 2 B'D + ca+A'B - a^a+AB' - ci-cBD' = 0, 

(10) b-b+C'D + c?C'A + ca+D'C - a^a+DC' - aSCA' = 0, 

(11) b-a+BD' + cb+DB' + b-cAB' - a^b+D'B = 0, 

(12) b-a+CA' + B+AC + b^cDC' - a^b+A'C = 0, 

(13) b-a+A'B + cb+B'A + b^cB'D - a^b+B'A = 0, 

(14) b-a+D'C + cb+C'D + b-cC'A - a-b+C'D = 0, 

(15) a-b + [B,B'] = 0, 

(16) a-b+[C,C] =0, (27) 

where we used the notations a_ = a(u/v), a + = a(uv) and similarly for b. Also A — A{u), A' = A(v) and similarly 
for B, C and D. Instead of writing these equations for K + , we solve the equation (9). The matrix M is determined 
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as the unit matrix 

M=(J J). (28) 

Therefore solutions for the boundary operator K + of the dual reflection equation is easily deriven. It is interesing 
to note that these sixteen algebraic equations for K_{u) do not depend on the twisted parameter t, i.e. there is no 
difference from the untwisted case to the twisted one related to U^s^)- Therefore the class of solution for Eq.(27) 
belongs to one analysed by Baseilhac and Shigechi [40]. 

The most simplest solution for the above equations consists of c- number elements, which is given by 

Kl{u)=( T_ + , + 2 U ~\ ^"-T 2) V (29) 
\ — [u — u ) ue_+u e + I 

where e± and k± are free parameters of the theory [13, 51]. Similarily, c-numbcr solution of the dual reflection 
equation is 



K +M~ { K-iq 1 ' 2 + q-V^qu 2 - q^u- 2 ) qV 2 UK* + q- x l 2 u~ X K J' (30) 

which parameters n, k* and n± are also free parameters. 

4.1 N — 1 dressed solution and Askey- Wilson relation 

To construct general TV-dressing i^-matrix, let us start from deriving the N = 1 dressed solution for the reflection 
algebra. By using the operator L(u) instead of the inverse operator of L(u), the N = 1 dressed solution 

K^(u) = LxHOJTIWMu-V) = ( ) (31) 

is obtained as the following forms 

A™ (u) = ue^ + u" V 1} + (u 2 - u- 2 ) {qWuW^ - q-^u^W^) , 
- ue W + u-h^ + (u 2 - u- 2 ) (q l ' 2 uW[ l) - q-^u^W^) , 



(1) (u 2 -vr 2 ) f k + k- W { yw^(q^ 2 u 2 + q-^ 2 u- 2 ) , G< 1J 

B (u) — 



k-w 02 \ c q 1 + q 1 



r w M (" 2 -^ 2 ) ( k+k^w^^/^ + q-^u- 2 ) G^ (1) » 

W" k+w W { c +,i/ 2 + 5 -i/2+^o J, (32) 

where and tu^ represents the Casimir operators (23) of the extended Sklyanin algebra (22) for j-th Lax 
operators. The parameters and uj^ are given by 

(1) _ k+k.w^iw^q-^vf + w^q^Vi 2 ) + Pe+e-W™ w 1 ? 

q-q 1 

Similarily in [27, 28], the generators g[ 1} and G^ are given by g[ 1} = [W^ , W (1) ], and g[ 1] = [W^ , w[\ in 
terms of the g-commutator 

[X, Y] q = q 1 ' 2 XY - q-^ 2 YX. (33) 
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Explicit representations of the generators w[ X \g^\g^ are written by 



TJ7 (1) 1,1/2, -1- ~ + , 1-1/2, . — | 



C 



(1) -4 2 V(g 1/2 +g- 1/2 )fc 2 -2 4^(g^+g^)M- , _ 1/2 2 ._ 1/2 - 2r+ , 2 ^- a) a) 

(l), ,-1/2, i/2 . _i/2w 1/2 -l - ~+ . -1/2 - — x . k + k^w^(q- 1 / 2 w ( L ) v 2 1 +q 1 / 2 w ( ^ ) v^ 2 ) 
-Wojk-h (g +g )(-g u i e-n2rf+g 7 Uie+r^Ti ) + (1) 1 , 

G i = z r 21 + (q ' Vi (t 2 ) +q 1 v 1 (r 2 + ) ) + cw K 0l 'w K Je-e+ 

k+k-wW (q-^w^v 2 + qVlw^Vi 2 ) 



-w, 



(1) Mi /2 (g 1/2 + g- 1/2 )(-g 1/2 «r 1 e+75if + + g-VVe-Tbifa") + 



■ 0l «-+«-! vy ti( a - y "a c +'2i< 2 ~r y fic-72i/ 2 ; "i rn — m ' 

Straightforward calculations show that the fundamental generators and satisfy the Askey- Wilson 

relations [52] 

[wf \ [Wf \ W (1) ],],-i = (g 1 / 2 + g" 1 / 2 ) 2 ^ 

+(g - g-^V^ - (g 1/2 + q-^k+k.w^w^e^, 

+ (q - q-^W^ - (q 1 / 2 + q-^k+k.W^ , (34) 

which lead to the g-Dolan Grady relations 

M^M^K^ [W^\[W^\[W^\w^%] q -,}^ P [M! ) 1 \wn (35) 

with pfp = = (q 1 / 2 + q^/^k-k+w^w^. For generic q, the generators are connected with the 

g-Racah polynomial and some related polynomials of the Askey scheme. 

The transfer matrix for the N = 1 deresscd solution can be expressed in terms of these four generators 

= JT« + (u 2 - u- 2 ){qu 2 - q- 2 u- 2 )l^\ (36) 

where 

4 1] = (g 1/2 + q- 1,2 W-e { l ] + K£ a) ) + (q^u 2 + q- 1 ' 2 ^ 2 )^ + K*e W ) 

-(q 1/2 + q- 1/2 )(u 2 - u- 2 )(qu 2 - g" 1 ^ 2 ) ( k + k - W j W + ] ( g l/2 u 2 + g-1/2,,-2) + w (l)\ / 



and 



^ = ( kWq + - r^iI) G i 7^71) G i I ■ ( 37 ) 

""o 2 



fc+w 0i fc_ 



The algebraic part I\ does not depend on the spectral parameter u, so that it is expilicitly separeated from the 
functional one. Then the eigenvalue problem of the transfer matrix is read in a different way as its problem of 
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X\. Realizing the algebra (22) by the difference operators, the eigenvalue problem of the conserved charge X\ 
leads to the second order difference equation. One can observe that eigenfunctions of X\ are associated with the 
g-hypergeometric function. The roots of polynomials, which obey the Bethe-Ansatz equation, is used to determine 
the spectrum of the system. We do not penetrate further in the details. Instead, we recommend you to read the 
references [53, 54]. 

4.2 iV-dressed solution 

To construct general iV-site lattice models, we would like to find explict representation of TV-dressed solution of the 
reflection equation. After the above calculation for the N = 1 dressing method, the same kind of calucuration 
is performed for the N = 2 case which result we suppress here. Based on these results and our previous results 
[27, 28], we can obatain the following form of the TV-dressed solution 

K^\u) = L(uv N )...L(uv 1 )Kl(u)L(u- 1 v 1 )-..L(u- 1 v N )=(^ ^("j )' ( 38 ) 

where these matrix elements are 

AW(u) = uef > + u~^_ N) + (u 2 u-2) ( uq^ J2 P -kW-k ~ ^ V 1/2 E P-MIl I - 



AT-1 JV-1 

fe=0 fe=0 
AT-1 N-l 



pW(tt) = ue {N) + u-h™ + (u 2 - u- 2 ) ( uq 1 / 2 £ P-V^+l - u^q^ 2 £ P^W (N) 



fc=0 fc=0 
N-l 



fe=0 



C (N \u) 



fc=l 
(u 2 - u- 2 ) 

N 



+ gl/2+ 1 g _ 1/2 E^ ) WC ) 1 ) • (39) 



fe=i 



Here the non algebraic parts -Pi^P, e^±\j^ N \u) and are given by the following equations: 



i N_1 /„l/2„,2 I „-l/2„,-2 \ ™ fc 

p(Af) = 1 I q ' u +q i u \ {N) 

-k " _l/2 i „-l/2 I „l/2i„-l/2 / " 



g l/2 + g -l/2 I g l/2 + g -l/ 

fc + fc_n(-4f)(-^) 

^V) = fc =) a _ tf'V + q-^u- 2 ) P r{u) + J Q N) 

AT 

(gV2_g-l/2) 



4 W) = (-^T^Fv^ ft «*(-»S ) )(-»£ ) ). 

fc=l 



with 



AT 

.(AT) 



C^ = (-l) JV -V /a + «" 1/a ) E 

fel<---<feiv_n-l=l 
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q -l/2 w W v 2 + q l/2 w W v ~2 c (°) c (0) ( 1/2 _ -1/2)2 



(« 1 /2 + g -l/2 )u ,g) lu g) 

Reccursive representation of the algebraic parts W^+i , Wl^P , , is given by appendix A. The form of 

the boundary K matrix (38) is directly proved by mathematical induction dressing from the N — 1 dressed solution 
to the N dressed solution. Then as the consitency condition we find the generalized linear combinations 

( 1/2 -\I2\ ( N ) N 

- {q N ~ q K — wiv + E ^ZMV-t + = 

fc + fc_n(-^ ) )(-4" ) ) 

fc=l 

(g V2 - g' 1/2 )4 jV) Ty(Af) i y r (N) W (N) (N) _ n 



fc=i 

(g 1 / 2 - q- 1 ' 2 )^ (at) A (JV) (A T) _ 

k=l 

{q 1 ' 2 - q- 1/2 )J N) - (JV) A (jV) ~ (A r) _ 
a7 '+1 + ^-Hl^HI+l _ U ' 



k=l 



The above relations for the operators are a natural consequency in order to form closed algebra in finite dimensional 
case. The transfer matrix is drived as the following form: 

tW(«) = Tr {K c + {u)K {N] \u)) 



N-l 

= T^H^ + ^-u-^-q-^^P^M^ (40) 

fe=0 

where the function (u) is 

J* N \u) = {q 1 ' 2 + q- 1/2 )( K *e { + N) + ne™) + (g l7 V + f^V 2 )^ + «*e^ } ) 

JV AT 

+ ( 9 V2 + ,-V2 )(u 2 _ u -2 ){qu 2 q -l u -2^ K+/ ^ + JJ(_40)) + n(-<)))^ (Ar) (-). (41) 



(=1 



and the algebraic parts I^fc+i are derived as 



By the construction for the transfer matrix, I^k+i snom d commute with each other. It turns out that the commu- 
tativity of l^k+i ^ s ensure d by the following algebraic relations 

[ W [fMV]=0, =0, ICC^O- [CIO=0 (43) 
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[w[ N k \ <?] = [<»<»], [<\ Gffl] = Gffi], G™] - [W^5, Gffi], (44) 

[Wffl. - Wi, Gffi], [Wffl, Cl = Wi. = IG&i, Gffi], (45) 

with 

[wf ) _ = ^ _ W W [14 ,(-) _ W W,G<^ fl _ l[Wi W - <»Gf] rl , (46) 

w° - c ^ - (gl/2 + f 1/2)3 ( n k* . ^ + • (48) 



Vm-l 



In fact these algebraic relations are verified by the mathematical induction with the reccursive expressions given 
in the appendix A. 

One can observe that the above algebraic equations possess the q-Doaln Grady relations even in the general 
N case. Some straightfoward calucurations show that G[ N ^ — [Wi,Wo] q and G[ N ^ — [Wo,Wi] q . From the initial 
conditions for the operator in the appedix A, we find the relations in the lowest order operators 



Pb 
1 

Pi 



w[v=wr+-^[wr,G[ N \, (49) 



where p^ = p[ N ^ = (q 1 ^ 2 + q l / 2 ) 2 k-k + Y\k=i w - ^ Therefore one can easily obtain the g-Dolan Grady 
relations 

[w[ N \ [w[ N \ [wrM N \] q -A = P ^[w[ N) M N \ 

from the algebraic relations in Eq. (43). 

5 Generalized McCoy- Wu model with open boundary conditions 

Here the generalized McCoy- Wu model with general open boundary conditions is considered as the most simplest 
model of our construction. The Lax operators L(u) and L(u) are expressed by 

L \ U )- T a *l/2=_ „l/4„-a,/4„.-l „-1/4„<t,/4„, i V> Z ) 



t 1 / 2 ~ca+ q^q-^u- 1 - q-^q^u 

through the fundamental spin-i representation on a symmetric relization 



& = Tq Tl/4 q T ° 3/ \ ^ - T9 T1/ V CT3/4 , r 12 = C( r_, f 21 = ca + , (53) 



where a± = (<7i±z<72)/2, and a.- tl (i = 1, 2, 3) are the Pauli matrices. Then the Casimir operators (23) is determined 

as 

w±=q Tl/2 , w 0l =w 02 =-l, w = (q + q~ 1 ). (54) 
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In addition, we set = l(i — 1, • • • N) to consider homogeneous spin chain model. The Hamiltonian Huw of the 
system is obtained from dlnt< rftt , where the transfer matrix is 



t (w )(u) = Xro (kc +{u)l (mw ){u) . .. L (w )(u) ^c (u )£(^) (u ) . ..Z#™0( U )) . (55) 
Differentiating lnt(u) with respect to the spectral parameter at u = 1 



du U ->1 yiq 1 / 2 + q- 1 / 2 ) (gV2- g -l/2) ^ ( g l/2 _ g -l/2) 

we can obtain Hamiltonian for the generalized McCoy- Wu model with general open boundaies as follows: 



*;=i 



+ (e+ + e -) l^T-^ + (g i/»-^ ^+ + J ■ (57) 

where A = (q 1 / 2 + q~ 1 / 2 )/2. The detailed analysis of the above derivation is given in the appendix B. It is clear 
that the Hamiltonian (57) for the untwisted limit, i.e. ti = 1 (i = 1 • • --/V), reduces to the XXZ open spin chain 
with general boundary conditions. By differentiating the transfer matrix n-times, mutually commuting quantities, 
which commute with the Hamiltonian, can be identified, although it seems that the method is much complicated. 

Alternatively, the transfer matrix by applying the Sklyanin's dressing method is described in terms of generators 
of the g-Onsager symmetry by taking the spin 1/2 representation (53) and its related Casimir operators (54). We 
denote the conserved quantities I^k+i m ^ nc s P m representation as ^fc+i- The Hamiltonian (57) is expressed 
by a linear combination of conserved quantities I^k+i fr° m our construction for the transfer matrix. Thus all 
integrals of motion in this model are explicitly determined as hk+i, i-e. one can show that 

[HmwJSIi} = 0, (k = 0, 1, • • • N - 1). (58) 

In particular, the conserved quantity l[ N ^ is given by 

= [nW { N) + n*w[ N) + ^-G[ N) + ^Gf >) , (59) 

where 

W[ N) = (fc+#V + fc_^ 1/2 a_) ® + ® Vt^ 1 '. (60) 

with the relations for g[ N ^ — [Wi, Wo] q and G[ N ^ — [Wo, W\] q . Instead of cosidering the eigenvalue problem Hmw, 
it may be read as the eigenvalue problem of the conserved quantity x[ N \ 



6 Discussion 

We consider the ancestor model related to the twisted V q (sl2) integrable models with integrable boundary condi- 
tions, which symmetric and nonsymmetric realizations of the algebra can generate several descendant lattice models 
without limiting procedures. The reflection equation and dual reflection equation for K±{u) docs not depend on 
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the twisted parameter t(ov 9). As expected that the same current algebra analyzed by Baseilhac and Shigechi 
is derived, we observed that the transfer matrix of the model is generated by the g-Onsager symmetry, and we 
identified all of fundamental generators of the g-Onsager algebra in terms of the extended trigonometric Sklyanin 
algebra. Although there is no difference between the untwisted model and the twisted model in the reflection 
and dual reflection equations, their generators have explicit dependence on the twisted parameter. The twisted 
parameter dependence for the generatros shows that the new quantum integrable model with integrable boundary 
conditions were explicitly constructed. 

We introduced the generalized McCoy- Wu model with general boundary condtions. Mutually comutting quan- 
tities of the genelized McCoy- Wu spin chain model with general open boundaries were explicitly expressed in terms 
of abelian subalgebra of the q-Onsager algebra. Therefore we confirmed that the model enjoys the g-Onsager sym- 
metry. To solve the spectrum problem is generally difficult in order that algebraic Bethe ansatz techniques have 
failed, because there is no reference state in the case of the general boundary conditions. Nice structure of the 
g-Onsager algebra may give useful information to solve this problem in similar to XXZ spin chain [55]. However 
this analysis goes beyond the scope of this article. We plan to do this problem elsewhere. 
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7 Appendix A 

We give the recursive representation for the operators , , and in terms of the untiwsted 

extended Sklyanin algebra by dressing the N — 1 dressed K matrix again. The direct method to find out TV-dressed 
solution for boundary K-matrix finds that 



(N) i -1/2 (AO 2 1/2 (AT) -2\ 



JV-l 

m_(«v2+ g -i/2)2 n(-<°x-4?) 



t-i N ~ 1 s 

' • i/2„-i- .,Sk)ssMN-i) 



-A_^ 1/ Vn2T 2 + ®(II(-< ) ))^ 



fc=l 



will - * o^K> - <r 1} ) + ^tt^ ® wit? - ^ ll2 ^:^f W) i * wr* 



fc + A_(«V2 + g -l/2)2 JJ(-4«)(-«,g)) V 



,_, JV-l N 

1/2 -i~ ~+ ^ i TT ( _n,S k )\\n( N ~ 1 ) 



-k + 4 V^if+ ® ( n {-<>))Gt 



k=l 
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N 

k + (qV 2 + q-V 2 ) H(-w£>) 



fe=i 

AT 



+^V /2 + <r 1/2 ) n(-o" ) ) (^'VW ® (w<K> - <r 1} ) - ?- 1/2 ^n 2 fr ^ (w f - 
fe=i 

V /2 + <T 1/2 ) 2 



and 



JV 



k + t N ]J( W 0l ) n l/2 v -2 ( ~ + ]2] 

^fe+l - — T 21 + / 1/2 +g -l/2) ^ ^ + ^Oi W 2 J ^ ^fe+l 

fc+(g l/2 + g -l/ 2) 

fe=l 

JV 

+MW 72 + ,-v2 } n( _^) } ( g vv w ® - 1J - - <r 1/2 ™r 2 - - <r 1) )) 
fc=i 

with the initial conditions 

AT 



fe=0 



fc + fc_( (Z 1 /2 +(Z -l/2 ) 2-Q ( _ u7 W )( _ w (^ 



fe=0 



fe=0 (gV2-g-l/2) 



and 

w(°) - c 

(-)" "z+i — c (-)i+i ' ^j+i — ^;+i 

8 Appendix B 



W*?=e<°> <>=e<°> GfJ^GffisefeW^-g-V), (i = 0, 1,2, • • •). (62) 



In the Appendix B, the Hamiltonian Hmw of the generalized McCoy- Wu model [41, 44] with general open bound- 
aries is explicitly deriven, because it is some complications compared with one of the XXZ opens spin chain with 
general boundaries. We can perform to differentiate \at{u) with respect to the spectral parameter u at u = 1 by 
using the Libnitz rule: 

d\nt^ MW \u) 

U—>1 

Tr (Kl(l)L^ w \l) ■ ■ • L[ MW \l)Kt(l)L[ MW \l) ■ ■ ■ lf W \l)) 
+Tr (Kl{l)tf W \l) ■ ■ ■ L[ MW \l)K c _(l)L[ MW \l) ■ ■ ■ L^ W \l)) 

+ £ Tr (K% (1)4 W) (!)••• L\ MW) (1) • • • L[ MW) (l)K<L{l)l[ MW) (!)■■■ L^ IW) (1)) 



du 

1 



£(MW0(1) 
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+Tr ( Kl(l)L { £ IW) (!)••■ L[ MW) (l)kl(l)L[ MW) (!)••■ L^ W) (1) 



AT-l 



+ ^ Tr ( tf°(l)L<f w > (1) • • • L\ MW) (l)K1(l)L\ MW) (1) • • • Z; (1) • • • L^ IW) (1) 



;(MW), 



i(MW) 



+Tr (1) • • • LP " V ' (l)K c _(l)Lr VV, {l)---L N ' (1) 

It is convinicnt for the derivetion that L( MVl/ )(l) and Z/ MW/ )(1) are expressed as 



- (MW) 



(MW) , 



x. (MW) 



- (MW) 



(1) 



+*; 



-1/2 







1/2 1 + ^ + 



(0 



-(MW), 



Vn tW 



/ 1 + 4° , ,1/2 1 " ° 



£(W) (1) = g(T W ) -lp 



(0 

3 



V 



o 

1/2 1 + ^ 1 ~ Q-3 ' 

tl 2 + 2 



where 7^ j denotes the permutation operator with properties 



for any operator Fj^ acting on lattice sites j and k. The first term is 



^^J^Tro (kUWT (I) - 4 W) (1)^- (D^ W) (I) - 4""° (1)) 

(gl/2 + g -l/2)( K + K *) r H + U J gl/2 + g -l/2' 



with the relation 



f< w )(l) = 2^(^/2 + g -V2)( c+ + e _)( K + K .) j. 
The second term contributes the left boundary term of Hmw- 

1— Tr (^(1)^(1) • • • M MW) (1)^(D^ W) (1) • • • 4 W) (D) 



1 



2(q - q- 1 )^ + k*) 
1 



Tro ^(1) 



2 AN) 



o 



t-2 2 W 



f -l/2 W 



1-cr 



(JV) 



2c(k + K*) 
A 



Tr n 



g V2 K + g -l/2 K . K+(<z l/2 + ? -l/2 )2g 



1 + a 



(N) 



4/; 



AN) 



7 l/2_i_„-l/2 



4i 



1/2 



1 - cr 



1 



K — K 



l/2 +g -l/2 (T + 



(JV) 



(63) 



(64) 



(gV2_ g -l/2) ( K + K *) 

where we used notations £ = (g 1 / 4 + g -1 / 4 ), r? = ((J 1 / 4 — g~ 1//4 ). Similarily, we find that the last term is equal to 
the second one. To obtain the third term, a key point relation is the followings: 



L ( MW) {1)L (MW) {l)i (MW) MW) (1) 



c (mw) (u ( £ 2 +V 2 vi l) 



t l/2 (0 



tr 1 /2 (7 W 
i- CT <*' 
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g3 (g l/2 + g -l/ 2) / l±^_ + h±l 



^ t-Y 2 ( i-^ +1 ') 



u 2^2 

l + \ ( l±^l + Qi-^) 

•3 q 1 /2+q-l/2 - „ 2 ' 2 U 



g l/2 +9 -l/2^+ "3 
-l/2 a _„(i+l) 

"2 



c -3 ( ,i/ 2+5 -i/ 2) / i±4^ + ^/ 2 ( i -^ +1) ) 



o , 3 j + 



2 +2 
V2 ,,1/2 



U 2^2 



Therefore the third term leads to the bulk interaction term: 

N-\ 



(0 toj+ftV 2 (i+i) m\ f 1 

+ g l/2 +g -l/2 ff - H 



(1/2 _ -1/2) 1/2 _ 

i= 1 

and the fifth term gives the same contribution as the third one. By using the relation 
L[ MW \l)k c _(l)L[ MW \l) 



(!) t-V2 n (1)^ 

\ _ / e+ - e_ 4fc+/c 



^ I 2 2 f V2„ J .JlK 1^01 



" 1 «; /2 d+4 l) ) , i-^ 1 ' i ri)j v ^ /'"■ 



XPoi 2+2 



( -l/2 f1 ,„(!)-, i 

Q *l (1+ CT 3 ) _|_ l-°3 



1+4" . t- 1 /^!-^ 1 ') 



2 2 

n «r 1/2 (i+^') , i-^ 

u 2^2 



= s 2 ((e + - + fti'V* 1 ' + *r 1/a *-^)) (J ?) 

the middle term can be identified with the right boundary term: 



1 



t( Mff )(l) 



Tr [K c + {l)L^ W \l) ■ ■ .L[ MW \l)K c _(l)L[ MW \l) ■ ■■L^ W \l) 



Gathering all of these results permits us to derive the form of Eq.(57). 
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